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Abstract 

A usual step in solving totally Schrodinger equation is to try first the case when di- 
mensionless position independent variable w is large. In this case the Harmonic Oscillator 
equation takes the form ( d 2 /dw 2 — w 2 )F = 0, and following W.h .B. method, it gives the in- 
termediate corresponding solution F = exp(-tu 2 /2), which actually satisfies exactly another 
equation, (d 2 /dw 2 + 1 — w 2 )F = 0. 

We apply a different method, useful in anharmonic oscillator equations, similar to that of 
Ram pal and Datta [1], and although it is slightly more complicated however it is also more 
general and systematic. 


After some arrangements Schrodinger equation for a simple harmonic oscillator is set as in (1): 


(d 2 / du 2 + c-o - C 2 U 2 )N{u) — 0 


( 1 ) 


h 2 c.Q = 2m E h 2 c 2 — {mu; 0 ) 2 (2) 

where c 0 and c 2 are the parameters of the differential equation, u — x — a is the distance from 
the particle to the point where potential energy is a minimum, u 0 is the classical angular frequency 
of the movement, E is the total energy and N is the probability amplitude for the particle to be 
found at u . In a very common procedure, (I) is first transformed to equation (3), in which w is 
a dimensionless independent variable, obtained by the mathematical manipulation (4): 

( d 2 /dw 2 + b - iv 2 )f(iv) = o (3) 


w 2 h = mwoii 2 hu>ob — 2 E (4) 

Some authors [2] get the solution of (3) for large id as a decaying exponential function of m/2, 
by means of the W.K.B. method, which is a factor in the total wave function. 

But others [3] simply propose the change ol dependent variable (5) without mentioning the 
W.K.B. method: 
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/ = exp{—w 2 /2)P(iv) (5) 

And a. few ones [4] consider that the dominant term, when w is large, is w 2 and they write (6) 
with b — 0 as an approximation of (3); and assert that (7) are its solutions: 

(d 2 /dw 2 — w 2 )F = 0 (6) 

F — exp(±tu 2 /2) (7) 

But after choosing the negative sign in (7), for the good behaviour of F , it is easy to demon- 
strate that (7) does not satisfy (6), and rather satisfies (3) when 6=1. We will use a method 
whose intermediate steps at all stages are correct. 

Starting with (1), the Ansatz (8) is composed of two factors: P that gives information of the 
zeros of N . and G: 


N = exp(-G)P (8) 

Non linear equation (9) is obtained from (8) and (1), with unknowns P and G : 

[P" - 2P'G') + \{G') 2 - G"]P = [-Co 4- c 2 u 2 ]P (9) 

[(G^_G'”] = -co(0) + c 2 ^ (10) 

( ’ means derivative with respect to a ). 

If P is an n degree polynomial P n , then for n = 0, Pq is a constant, and (10) is a non-linear 
equation with only one unknown. By watching (10) it is noticed that (11) is the solution of (10) 
if constraints (12) hold: 


G = (3u 2 (11) 

4/i 2 = c 2 20 = co(0) (12) 

The first eigenvalue £(0) can be obtained from (2) and (12). Consequently, (13) is the solution 
of the Schrodinger equation for n = 0; 0u 2 as the argument of the exponential function must be 
dimensionless: 


N = / 0 = P 0 exp(-/3ti 2 ) 


(13) 


Therefore iv becomes a dimensionless variable and (1) is transformed into (7), and taking into 
account (12): 


w 2 = aftii 2 

(14) 

[d 2 jdw 1 -1- c Q (n)/a0 — 4 w 2 /a 2 \f n {w) = 0 

(15) 
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(15) is precisely (3), with b{n) = c 0 (n)/a/i, and 6(0) - Co(o)/afJ = 2/a. (16) is the correct 
solution of (15) and (17) is the differential equation of the polynomial P n : 

N n = P„exp(-u; 2 /2) ( 16 ) 

[rf 2 /^ 2 ” (4wfa)dldw + 2(6 - 1 )/q]P„ - 0 (17) 

For a = 2 and 6(?i) - 2u + 1, (17) is the equation of Hermite, and P n = 7/ 7l are Hernnte’s 
polynomials. The energy eigenvalues can easily be obtained. 

The Simple Harmonic Oscillator Schrodinger equation is perhaps the known differential equa 
tion with the most accurate solution. It would not be worthy to obtain that solution again, if 
methodological aspects are not taken into account. If a physicist plans to work in problems related 
to differential equations, it is useful to give her ( him ) general and powerful methods. As the 
Simple Harmonic Oscillator is the first approximation to many physical models, and one of the 
first problems with which the students are put into contact, it is good to take advantage of meth- 
ods that can be used in better approximations to more complex physical models, and more exact 
formulations, as relativistic ones lor instance. The author has made a review of the relativistic and 
non- relativistic isotropic harmonic oscillators, and uniform magnetic fields [5], using this method. 
With polar coordinates and centrifugal potentials, other polynomials depending on two quantum 
numbers, and one extra factor are the solutions of the radial equations, first the author had 
used the method in solving anharmonic rectilinear oscillators, and anharmonic isotropic oscillator 
equations [6], and continues working further these topics. In all those cases mentioned, and heie, 
we consider that writing equations (9) and (10), is the most important step that permits to find 
the solution for large uk and the independent dimensionless variable. 
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